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The density-functional (DF) theory provides a simple method for calculating the proper- 
ties of an interacting system under an external potential by associating it with a correspond- 
ing non-interacting system. Here, we find some relations in this non-interacting system, 
which can enable us to set up the thermodynamical relations for a neutral electron-nucleus 
mixture in terms of quantities of the non-interacting system and the exchange-correlation 
effect. In this way, Andersen's force theorems are clearly described and easily proved in con- 
junction with other types of force theorems, Janak's theorem extended to a finite-temperature 
system, pressure formulas and some thermodynamical relations. For this purpose, thermo- 
dynamics in the DF theory is presented in a systematic and explicit way. 

§1. Introduction 

The density-functional (DF) theory has been used with success to investigate 
the equilibrium properties of solids and liquid metals in the almost cases when the 
electrons can be treated as being in a perfectly degenerate state. Therefore, not so 
much attention is paid to thermodynamics in the DF theory, although the DF theory 
is extended to be applicable to finite-temperature systems As a consequence, 
we found some confusions in applications of the DF theory to condensed matter, 
as shown by some examples below. Frequently,^ the pressure of a solid is 
determined by the derivative of total electronic energy E$ with respect to the atomic 
volume Q (a uniform compression of the lattice): P = —dE^/dQ. To be correct, 
the volume derivative of electronic energy should be done under the fixed nuclear 
positions (fixed lattice) for obtaining the 'electron' pressure, as shown in the present 
work. For a next example, the proofs of a generalized Janak's theorem to a finite- 
temperature and the familiar Fermi distribution law are inappropriately proceeded,^ 
since the thermodynamic relations in the DF theory are not recognized clearly. The 
last example is the entropy expression^ in the DF theory, the derivation of which is 
not transparent. In this situation, it is useful to write up thermodynamics in the DF 
theory applied to a neutral electron-nucleus mixture in the explicit and systematic 
form. 

Andersen^ has derived force theorems^ <EH>E30 associating a certain change in 
the sum of eigenvalues determined by the DF theory due to the volume change 
or atomic displacements with the pressure or the forces, respectively. However, the 
expression of force theorems is ambiguous and not definite to be used for calculations 
of the pressure and the forces. In addition, their 
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basis of the physical picture as HeineP^ has done even though in long four pages. 
Furthermore, the force theorems are proved only at zero temperature and in the local- 
density approximation (LDA). Contrary to this situation, thermodynamics in the DF 
theory is shown to provide clear representations of the force theorems extended for a 
finite-temperature system without use of the LDA, in conjunction with their simple 
proofs in the present work. 



§2. Thermodynamics of inhomogeneous system 
under the external potential U(r) 

In this section, we sketch a brief outline of thermodynamics of the inhomoge- 
neous system for the purpose to set up thermodynamics in the DF theory. When 
the thermodynamic potential fl of an inhomogeneous system under the external po- 
tential U(r) is given, the density distribution n(r) is determined by the functional 
derivative of f2 with respect to 7(r) = \i — U (r)Pi 



<h(r) 



5tt 



TV 



5U(r) 



-n(r) 



(2-1) 



TV/i 



at temperature T of the system with a volume V and chemical potential /i. Therefore, 
the natural variables of the thermodynamic potential Q for this inhomogeneous 
system are T, V, j(r) or T, V, fi and U(r), and the change in Q associated with 
any infinitesimal variation in these natural variables is expressed by 



(2-2) 
(2-3) 



dQ = -SdT - PdV - J n{r)5 1 {r)dr 

= -SdT - PdV - Ndfi + J n(r)5U(r)di 



The relation (|2-lj) enables us to replace an independent variable 7(r) by n(r) with 
use of the Legendre transformation of il, which introduces the intrinsic free energy 
F: 

F = n-[J^-r j(r)dr = Q + / n(r)-y(r)dr . (2-4) 

TV J 



#7(r) 



Now, the natural variables of F become T, V, n(r), or T, V, N, <r(r), and the 
variation is given by 



dF = -SdT - PdV + J 7 (r)5n(r)dr 



(2-5) 



-SdT - PdV + 



H-j U(r)a(r)dr dN - nJ U(r)5a(r)dr . (2-6) 



Here, cr(r) is the shape factor defined by the relation n(r) = iV<r(r)!^ 

In a similar manner, the free energy F, the internal energy E, the intrinsic in- 
ternal energy E and the intrinsic thermodynamic potential Q for the inhomogeneous 
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Table I. State functions and their changes with respect to infinitesimal variations in the natural 
variables. 



state function 


natural variables 


thermodynamic potential Q 


T, V, 7(1-) or 


T, V, (i, U(r) 


dQ = -SdT - PdV - f n(r)<5 7 (r)dr 




(Al) 


= -SdT - PdV - Ndji + / n{r)5U{r)dr 


(A2) 


intrinsic thermodynamic potential Q 


T, V, n, n(r) 


dQ = -SdT - PdV - Ndfi - J U{r)5n(r)dr 


(B) 


free energy F 


T,V,N,U(r) 


dF = -SdT - PdV + fidN + / n(r)SU(r)dr 


(C) 


intrinsic free energy F 


T, V, n(r) or 


T,V,N,a(r) 


dF = -SdT - PdV + J 7(r)£n(r)dr 




(Dl) 


= -SdT-PdV+ [n-jU(r)a(r)dr]dN-NjU(r)5a 


(r)dr (D2) 


internal energy E 


S,V,N,U(r) 


dE = TdS - PdV + fidN + / n(r)5U(r)dr 


(E) 


intrinsic internal energy E 


S, V, n(r) or 


S,V,N,a(r) 


dE = TdS - PdV + / 7 (r)5n(r)dr 




(Fl) 


= TdS -PdV + [n-JU(r)a(r)dr]dN 


-NjU(r)Sa(i 


•)dr (F2) 



system are introduced by means of the following Legendre transforms: 

= Q + fiN , 



^ „ dQ 
F = Q-n — 



E = F-T 



E = F-T 



dF 
df 

dF 



TVU(r) 



VNU(r) 



F + TS. 



dT 



F + TS. 



Q = Q 



Vn(r) 

5Q 



6U(t) 



U{r)dr = Q- n(r)Z7(r)dr, 



(2-7) 
(2-8) 

(2-9) 

(2-10) 



of Q, F, F and Q, respectively. Note that 'intrinsic' thermodynamical quantities Q 
(such as E, F and Q) are related to the usual quantities Q by the relation: 



Q = Q+ U{r)n{v)dv. 



(2-11) 



Changes of these thermodynamic quantities associated with infinitesimal variations 
in the natural variables are summarized in Table. [|]with additional new relations. 
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Equations (|2-2|) . (|2-3[) . (C) and (E) of Table. [Ulead to formulas to determine the 
density distribution n(r) in several ways as 



5f2 



rar 



<57(r) 



TV 



5U(r) 



5F 



6U(r) 



5E 



TVN 



SU(r) 



(2-12) 



SVN 



Inversely, the external potential U(r) is determined as 



5F 



5n(r) 



5E 



TV 



5n(r) 



(jl-U(t), 



(2-13) 



sv 



by using the relations, (Dl) and (Fl) of Table [IJ Similarly, the chemical potential fj, 
and the pressure P of the inhomogeneous system are calculated by the formulas 



dF 



TVU(r) 



dE 
dN 



5F 



SVU(r) 



5n(r) 



5E 



TVU(r) 



Sn(r) 



SVU(r) 



and 



dF 




dE 






dQ 




dV 


TNU(r) 


8V 


SNU(r) 




dV 


T 7 (r) 


dF 




dE 






dfi 








Tn(r) 


~ dV 


Sn(r) 


W 


Tfin(r) 



(2-14) 

(2-15) 
(2-16) 



Also, the equilibrium density distribution n(r) satisfies the following condition: 



5Q 



5n(r) 



5F 



TVfJ,U{r) 



5n(r) 



5E 



TVNU(r) 



5n(r) 



0. 



SVNU(r) 



(2-17) 



§3. DF theory of electrons under the external potential 
caused by the fixed nuclei 

Let us consider an electron-nucleus mixture consisting of A^ n nuclei and ./V = ZN n 
electrons with an atomic number Z. When the nuclei behave as classical particles, 
the nuclei are considered as generating the external potential U (r) for the electrons to 
form an inhomogeneous electron systemff^ this external potential for the electrons 
is the Coulomb potential produced by the nuclei fixed at the coordinates R a : 



U(t) 



Ze 2 



R„ 



(3-1) 



In treating this inhomogeneous system, an effective external potential f7 e ff( r ) can De 
defined so as to satisfy the condition that the true electron density n(r\U) under the 
external potential U (r) should be identical with the non-interacting electron density 
distribution n°(r\U e s) under the external potential U e fi (r): 



n°(r|£4ff) = n(r\U) 



(3-2) 
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Here, the non-interacting electron density distribution n°(r\U e g) is determined as 

>|[/ cff )^/,|<Mr)| 2 , (3-3) 



n 



in terms of the wave function fa (r) and eigenvalues t{ , which obey the wave equation 
for a single electron: 



-h 2 

2m 



V 2 + U eS (r) 



(3-4) 



and the Fermi distribution /j = l/{exp[/?(ej— /io)+l} with the chemical potential /io 
of a non-interacting electron gas. 

For the purpose to obtain an explicit expression for the external potential U e g (r), 
let us note the relation for intrinsic free energy of non- interacting electrons ft(= F s ): 



8ft 



5n°( 



5ft 



TV 



5n(r) 



Mo - U eS (r) = 7 e ff(r) , 



(3-5) 



TV 



which is ()2-13[) rewritten for the non-interacting system. Here, the intrinsic free 
energy ft\n ] of the non-interacting system is written in an explicit form: 

oo 

ft[n°} = --^2ln{l + eM(3(^-e l )}}- J n°(r)[/ eff (r)dr + (3-6) 

= T s [n°]-TS s [n }, (3-7) 

where T s and S s are the intrinsic internal energy E s and the entropy of the non- 
interacting system, respectively, defined by 



T s [n°) = E s = E s 



n°(r)C/ cfT (r)dr 
4>*(r)=^V 2 fa(r)dr. 



S s [n°] = -k B Y}h^fi + (1 - /i)ln(l - / f )] 



with the internal energy E s = ^ fiU- On use of the relation: 

8F 



5n(r) 



fi-U(r) =7(r), 



TV 



the effective external potential f7 e ff( r ) is represented explicitly as 

5Ji 



C/ eff (r) = C/(r) + 



8n(r) 



Ml 



(3-8) 
(3-9) 
(3-10) 

(3-11) 
(3-12) 



TV 



with ft = F — F s and /xi = /J — /io, which is derived by subtracting (|3-5[) from (j3-ll|) . 
Equation ()3-12[) is the effective external potential in the framework of the DF theory. 
At this point, we introduce the exchange-correlation free energy by 



(3-13) 
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in terms of the intrinsic electrostatic energy E es , a part of the electrostatic energy 
E es [n] described as 

r , e 2 f nfr)n(r') , . 1 ^ Z 2 e 2 A f Ze 2 n{r) , , . 



a^tp ^ 1=1 

= E cc + E im + .Ben = E cs + y n(r)E7 (r)dr . 

Then, we can write the intrinsic free energy F in the form: 

F = T s + J\ = T s + F xc + E cs , 



(3-15) 



(3-16) 



within the framework of the DF theory. 

With use of the Legendre transformations of F, we obtain the thermodynamic 
potential in the form: 



Q = F 



5F 



6n(r) 



n(r)dr = F — fiN 



where 



TV 



1 

Q s = -J2 ln(l - h) = — ln ^ + & XP\P^0 ~ *)]} , 

i i={ 



O — f 

J ^xc J xc 



5 Fx, 



Sn(r) 
SE„ 



n(r)dr . 



TV 



5n(r) 



n(r)dr = E nn - E ee , 



(3-17) 
(3-18) 

(3-19) 
(3-20) 
(3-21) 



TV 



and the intrinsic internal energy E in the form: 

OF 



E = F — T 



dT 



T s [n] + E xc [n] + E a 



(3-22) 



Vn(r) 



with the exchange-correlation part of the internal energy defined by 

dF xc 



E„. In] = T xc - T 



dT 



Vn(r) 



(3-23) 



From the above expressions for the intrinsic quantities, E and F, the relation (|2-11|) 
provides the internal energy E: 



E = E+ n(r)U(r)dr = T s + E xc + E es , 



(3-24) 
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and the free energy F: 

F = F + J n(r)U(r)dr = E-TS 

with the definition of the entropj^'^ 



TS S + JF XC + E c 



S = 



dF 
dT 



VNU(r) 



dF 
dT 



Vn(r) 



dT 



Vn(r) 



(3-25) 
(3-26) 

(3-27) 



Equation (|3-27|) is obtained from Eq. (j3- 16|) and the relation: dE cs /dT\ Vn ^ = 0; this 
derivation is different from those of the references,^" which are in a roundabout 
way. 

By use of the relation 



A*xc(r|nJ 



8n(r) 



5E^ r In] 



TV 



SV 



we can express the effective external potential (|3-12j) as 

6E es 



U eS (r) 



5n(r) 



U(r) 



5n(r) 
nj a; 

+ /u xc (r|n) - (ii . 



Thus, the effective external potential is rewritten in simple forms: 



U eS (r) 



6[E es + ^ xc ] 



5n(r) 



MI 



S[E es + E xc ] 



Sn(r) 



Mi 



SVU(r) 



(3-28) 



(3-29) 



(3-30) 



TVU{r) 

where we make a convention that the derivative 5E es /5n(r) is performed with the 
external potential U(r) (that is, the nuclear coordinates {R Q }) being fixed. 

§4. Thermodynamics of an ideal gas under the external potential [/ e ff (r) 

As shown in the previous section, the DF theory reduces the many-electron 
problem under the external potential to the one-electron problem by associating an 
interacting system to a non-interacting system (an ideal gas under the effective exter- 
nal potential). As a consequence of this association, the thermodynamical properties 
of an interacting system are described in terms of those of the ideal gas under the 
effective external potential. For this purpose, we enumerate here some important 
properties of an ideal Fermi gas. 

The change of the intrinsic internal energy E s = T s [n] of the non-interacting 
system is given by 

dE s = TdS s - P dV + J 7 eff (r)(5n(r)dr = dT s [n] , (4-1) 

while the change of the intrinsic free energy F s = T s is 

dF s = -S s dT - P dV + f 7 cfr (r)5n(r)cir = dF s [n] . (4-2) 
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Therefore, from the above equations, the external potential 7eff(r) = fXo — U c ff(r) is 
obtained by the functional derivatives in the two ways: 



5F a [n] 



5n(r) 



STJn] 



6n(r) 



Mo - U eS (r) = 7 e ff(r) , 



(4-3) 



s s v 



which leads to (|3-28|) with aids of (|2-13|) . Also, the pressure Pq is derived from (|4-1|) 
and P^jl as 



3P V = -3V 



v 



dV 



-3V 



n 



S s n(r) 



dV 



Tn{r) 



= 2T s [n] - J n(r)r • VU cS (r)dr 

= <f r • • dS , 

JdV 

where the kinetic pressure tensor is defined b^^^ 

V-P° K = -n(r)V*7 cff (r) . 



(4-4) 

(4-5) 
(4-6) 



(4-7) 



The above result is derived as follows. First, with use of the scaled density n\(r) 
A 3 7j(Ar)p3 the kinetic energy functional is shown to satisfy the relation: T s [n\] 
\ 2 T s [n], which yields 



-3V 



dT s [n] 



dV 



X 



dT s [n\] 



S 3 NU ef! (r) 



dX 



2T s [n] 



(4-8) 



A=l 



Next, we note that the volume derivative of T s at the fixed S s , N, U e g is expressed 
by a sum of the volume derivatives of T s at the fix n(r) and via the change in n(r): 



n 



dV 



dT s [n] 



S 3 NU eS (r) 



dV 



+ 



S s n(r) 



5TJn] 



Sn(r) 



dn[r) 



SsV 



dV 



dr , 



SsNU c g 



(see Appendix). Here, the first term in the right side of (|4-9j) is expressed as 



-3V 



dV 



STJn] 



S 3 n{r) 



2T s [n] + /n(r)r.V Mr) 

with help of ()4-8|> and the following equation?^ 

- / f(r)^dr = X^r f /(r)n A (r)dr 
J dV dX J v/X 3 



dr . 



(4-9) 



(4-10) 



s s v 



A=l 



A-4 / /(r/A)n(r)dr = - I n(r)r • V/(r)dr , (4-11) 
dX J v Jv 



which is applied to the second term in the right side of (|4-9[) . Finally, as a consequence 
of (|4-1U|) and (|4-3|) . we obtain (|4-5jl . which may be written in the surface integral 



Thermodynamics in Density- Functional Theory 



9 



1)4- 6j) over the surface dV of the system.^ It is important to recognize that the 
kinetic energy of a non-interacting electron gas is written as 

2T s [n] = [ n(r)r • VU cS (r)dr + I r • • dS = [ trP^dr , (4-12) 
Jv JdV Jv 

which yields the surface- integral expression ()4-6|) . The surface integral term is omit- 
ted in the standard treatment G3JIDJ3 which yields erroneous results if we treat a sys- 
tem with the non-zero electron pressure by using ()4-12j) without the surface integral 
term. For example, Slatei^ derived the virial theorem, 2(kinetic energy) + (potential 
energy)= — ^2 a ~R, a X7 a E, which is valid only for a solid with the zero electron- 
pressure, as will be shown afterward by ()6-8|) and (|6-1U|) . 

When we differentiate the internal energy E s = £^ fi€i in a simple manner, we 

find 

d £7^ = 6; + /^. (4-13) 



dfi 



dfi 



It is important to see the reason why the second term is discarded to yield in the 
circumstance that /j is a function of e^. As is seen later, this reasoning provides a 
clue to relating thermodynamics in the DF theory to usual thermodynamics. First, 
we consider the derivative of E s with respect to an arbitrary parameter A by noting 
the relation, E s = ^ /je, = Q s + TS S + fioN. Because of dQg/dT^^u^ = —S s , 
we get the following relation: 



Tx la ' + TS ° 



dfic 



d\ 
dfi s 



+ 



T 



dX 



+ T- 



T 



dT 

, dS s 
dX ' 



V/i U cl 



dT 



dX 



, dS s 
dX 



(4-14) 
(4-15) 



where the first term of (|4-15|) is put in the form 



dQc 



dX 



dQ q 



dX 



H 



TVU F 



d^o df2 s 
~dX~ 



dX 



Thus, from (J4-15|) and (|4-16j) we obtain a simple expression (see Appendix) 

dE s df} s 



dX 



dX 



dS s dN 
~dX ^°~dX 



If we take as A = /$, the above equation takes the form 



dE„ 



dfi 



vu cB 



+ W) 



which is consistent only with the following three equations: 

1 



TVlM>U eS 



ft = f{u- Mo) 



exp[/3(ej - no)] + 1 



(4-17) 



(4-18) 



(4-19) 
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d 

Wi 
dQ s 



vu cB 
= 0. 



(4-20) 



(4-21) 



"J * TV^ U eS 

By taking account of this situation (|4-17|l . we set up a formula to perform a 
derivative of intrinsic free energy E s = T s [n] by A, which is not contained in U e g, in 
the form: 

dT* df2 q 



dX 



dX 



This equation plays an important role to associate the DF theory with thermody- 
namics of an inhomogeneous electron system caused by the external potential of the 
fixed nuclei. The proof of this equation is shown as follows. The kinetic energy 
functional T s is written in the form 

T S = E S - j n(r)U cS (r)dr = f2 s +TS s + J n(r) 7cfr (r)dr , (4-23) 
and Eq. (|4-15|) is written as 



dX 



[n s + ts s 



dfi<. 



dX 



T 



, dS s 
dX 



2^ dX + 2^7x^-^ 



dX' 



Here, the second term of (|4-24|) is written as 



dfi. 



dX 



dQc 



dX 
df2 s 



+ 



dX 



Tfi U eS 



^7 eff (r) dn ( r ') drdr > 
TV 5n(r') dX 



^7eff(r) 

^7 cff (r) dn(r') 
5n{r') dX 



n[r) 



drdr' . 



Also, the third term of the right side of (|4-23|) generates 
d 



dX 



™(r)7 efr (r)(ir 



,dn(r) 

^ <ix~ 



57 cff (r) dn(r') , 

n(r) MrTir" ■ 



(4-24) 

(4-25) 
(4-26) 

(4-27) 



In final, a cancellation between the second terms in the right side of (|4-26|) and (|4-27|) 
results in (g22>. 

§5. Force theorems 

In this section, we derive three types of Andersen's force theorems on the basis 
of (|4-22p and (|5T6|) . The first type is a force theorem concerning the pressure. The 
pressure of a inhomogeneous electron system under the external potential caused by 
the fixed nuclei at {R Q } is defined by (|2T5I) . and can be calculated by introducing 
the scaled electron density distribution n\(x) = A 3 n(Ar) in the following form: 

dE v [n;{R a }] 



3P P y 



-3V 



dV 



X 



SN{R a } 



dEx 
dX 



(5-1) 



A=l 
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with 



E X = E y / xZ [n A ; {R Q }] = Tj^ [n x ] + [n x ; {R a }] + E^ [n x ] 

With aids of 032), we get^ 



dA 



SiV{R Q } 



dX 



+ 



T/loUeff 

dej 

dX t . 



V/A3 



n A (r)[^ - C4flf(r)]dr 



+ / n(r)r- VC/e ff (r)dr. 



Since the electrostatic energy (|3-15j) satisfies the following relation 

£ e 7 A3 [n A ; {R Q }] = A^[n; {AR a }] , 
the A-derivative of the electrostatic energy^ becomes 



A 



dEeJ X [n\] 



dX 



A=l 



-^es / Ra ' Fq. — / 



(5-2) 

(5-3) 
(5-4) 

(5-5) 
(5-6) 



where the force on a-nucleus is given by F c 
correlation energy 

p ar P» of (j52J) yields 



I \ Y7 ^ Ees A 

n(r)r • V - . , dr , 
on[r) 

V a E es . Also, the exchange- 



A 



dE^c X [nx] 



dX 



x=i 



trP x ,dr 



V 



n rjr 



v 



V— ^-dr +& r ■ P xc • dS , (5-7) 

<M r ) lav 



where the exchange-correlation pressure tensor P xc is defined by 



V ■ P xn = n(V)V 



8n(r) 



= n{r)V — - 

TV ° n ( T ) 



(5-8) 



sv 



Because of (|3-3U|) . the sum of the second term in (|5-4|) . the last term of (|5-6|) and 
the first term in the right side of (|5-7|) becomes zero owing to (|3-3UJ) ; as a result, a 
pressure formula (the force theorem concerning the electron pressure**) is written as 



3P e V = -3V^2f\ 



dej 
dV 



+ (b r • P xc • dS , 

u aB Jdv 



(5-9) 



in terms of eigenvalues e, of the Kohn-Sham equation l|3-4|) . 

Another important relation is derived by noting the relation^ 



TY^[nx] = X*Tj[n] 



which leads to 
dX 



SN{R a } 



2T.H = £ hji 



+ / n(r)r • VU cfi (r)dr 



(5-10) 



(5-11) 



/ trP^dr = <£ r • P^ • dS + / n(r)r • VC/ cfr (r)dr , (5-12) 
/ay Jv 
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with combined use of (|5-4[) . From the above equation, we relate the volume derivative 
of eigenvalues e, to the surface integral of the pressure tensor P° for a non-interacting 
electron gas: 

<] ' ' I r • P° K -dS = 3P V, (5-13) 

JdV 



dV 



which alters a pressure formula (j5-9|) in the form 



3P e V = <£ r-(P° K + P xc ) • dS . 
JdV 



(5-14) 



At this point, we remark that Eq. (|5T3|) is directly derived as the pressure of the 
ideal gas from (|4-6j) and the thermodynamic relation: 



3P V 



dV 



TfloU c a 



dej 
dV 



(5-15) 



The formula ()4-22j) to calculate the A-derivative of T s is valid only when the 
parameter A is not contained in the external potential U e R . When the parameter A 
is involved in the external potential U e g (for example, A = R a ), the formula (|4-22l) 
is changed in the form: 



dX 



dQc 



dX 



^dS s f , .dn(r) 

+ T df + / 7cff(r) ^T 

T W n(r) dA J dA 



dr 



n(r) 



dU eS (r) 



dX 



n(r) 



dr, (5-16) 



where the derivative of Q s is performed under the fixed density distribution n(r), 
instead of the external potential U e s being fixed as (|4-22j) . and the last term of (jfvlfij) 
contains the direct variation of U c r. Actually, Eq. I|5T6[) is derived in a similar way 
to get Eq. (|4-22|) with use of the fact that T s = fi s + TS S + fioN and that the natural 
variables of the intrinsic thermodynamic potential Q s are T, V,//o and n(r). When 
the A-derivative is performed under the fixed S s and V, Eq. ()5-16|) becomes 



dTs 
dX 



s s v 



H dX 



+ 



n(r) 



.dn(r) 



n(r) 



dU eS (r) 



dX 



dr . 



n(r) 



On the other hand, this derivative can be represented in another form: 

dn(r) 



dT s 




f ST s [n] 




dX 


SsV J 


' 5n(r) 


SsV 



dX 



-dr 



dn(r) 



(5-17) 



(5-18) 



owing to (|4-3|) , By comparing (|5-18|) with (|5-17|) . we obtain a formula to calculate 
the force on the atoms in the form: 



Ef ^± 
h dX 



n(r) 



ra(r) 



dU eS (r) 



dX 



dr . 



n(r) 



(5-19) 



When A = R a , this equation assumes the form: 



da 



n(r) 



n(r) 







dR, 



E 



-Ze 2 



■a I r -K-y 

7 " 



-Ze 2 



. j f ' ,' l "., r/r. (5-20) 
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Therefore, the force on a-atom is represented as 



d~R n 



n(r) 



(5-21) 
(5-22) 

(5-23) 



If we expand the eigenfunction fa by using a local base set itj(r) as fair) = J2j u j ( r ) c i ; 
the force is described by the matrix form as is given by 



dei 



dR r 



n{r) 



£(<* 



d(H - ei S) 
3R„ 



\Ci) , 



(5-24) 



with Hij = (ui\H\uj) and Sij = (■Uj|«j).™ 

Next, we consider another type of force theoremj*^'^ when not only the nucleus 
at R Q , but also the electronic charge within an arbitrary volume Q a enclosing the 
a-nucleus is displaced rigidly. Instead of displacing the volume Q a , it is more conve- 
nient to think in terms of the displacement of the rest system by — 5H a . Therefore, 
Eq. (j5-19j) assumes the form 



5ei 



n(r) 



n(r) — 3 -dr 



(5-25) 
(5-26) 



n(r)W cff (r)dr. 

On the other hand, the Maxwell tensor Tm obeys the following equation^ 

-V • (P° K + P xc - T M ) = Yl F « s (* ~ R -) • ( 5 ' 27 ) 

a 

Thus, integral of ()5-27(l over the volume fi a yields the force on a-nucleus in the form: 



/ (P°K + Pxc " T M ) • dS 

J dO a 

f V • (P° K + P xc )dr + 1 T M • dS 

n(r)VC/ e fj(r)dr — / n(r)V// xc (r)dr + <z> 
J Q a Jd 



Tm • dS 



d~€i 



<5R„ 



n(r) JO, 



ra(r)V// xc (r)dr + 



M 



dS 



dO a 



(5-28) 
(5-29) 
(5-30) 
(5-31) 



where the force is described in terms of eigenvalues e\ due to (|5-26|) . Equation (j5-31j) 
should be compared with (193) of the reference,^ where no exchange-correlation 
contribution is contained and the electrostatic force is differently defined. 
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§6. Thermodynamical relations in DF theory 



In previous section, we have proved Andersen's force theorems on the basis of 
(|4-22j) and (|5-16|) . In this section, we show that these equations, ()4-22(l and (|5-16|) . 
provide several thermodynamical relations in the DF theory. First, we show that the 
force F a on a-nucleus in the system is described by the thermodynamic functions 



-V Q E\ SVN 



-V Q F\tvn 



For example, the relation, F c 
following equations: 



tv, ■ (6-1) 
V a [T s + E es + E xc ]\svn, is proved by use of the 



V Q \E es + E- 



*c\\SVN 



Vc E e 



ln(r) 



5[E es + E x 



6n(r) 



SVU(r) 



V Q n(r)dr, (6-2) 



V Q t s \ svn 



5T S 



5n(r) 



V a n(r)dr 



s s v 



[^eff(r)+/ii]V Q n(r)dr, (6-3) 



which lead to F a = — V a E es \ n ^ with help of (|3-3Ujl . In the above, we used an 
important relation to cancel the term involving U e g in (J4-22|) and (|5- 16|) 



d[E cs + E xc ] d[E cs + E x 



dX 



r/A 



+ 



n(r) 



dU(r) 
dX 



n(r)dr + 



5[E es + E x 



5n(r) 



dn(r) 



SVU(r) 



dX 



dr , 



(6-4) 

which is derived by applying (|A-2}1 to the intrinsic internal energy E = E—J U(r)n(r)dr 
with the natural variables, S, V and n(r). At this point, it should be noticed that 
the relation (J6-3|) is equivalent to the following equation 



n(r) V a U eS (r)\ n ^r = V a f2 s 



du. 



n(r) 



r — R a 
n[r)- — — —dr = , 



R I 3 



(6-5) 
(6-6) 



which results from ()5-16(l . Therefore, the force expressions ()6-l(l are equivalent to 
(I6-6|) . that is, the force theorem ()5-23() . Similarly, other relations in l)6-ip can be 
proved with use of (|5-16|) , 

Equations (|2-15j) and ()2-16|) provide two different but equivalent expressions for 
the pressure as shown below. First equation (j2- 15f) generates a pressure formula in 
the form: 



TN{R a } 



5W{R a } 



2T S [n] + E CS + J trP xc( fr - ^ R 



(6-7) 
(6-8) 



Equation (|6-8|) is proved from (|5-6|) and (|5-7|) with additional use of ()5-ll(l . Here, it 
should be kept in mind that the electron pressure is determined by the derivative 
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with the nuclear positions (that is, the external potential U(r)) being fixed . In 
some pressure calculations p»J3Ja this fixed-nuclei condition is neglected; as a result, 
the nuclear virial term in Q6-8JI disappears in the pressure. On the other hand, the 
second equation (|2T6|) leads to a pressure formula in the form: 



sp e y 



-SV 



dF 
dV 



-SV 



Tn(r) 



dE 
W 



5n(r) 



dv 



r • (Pk + Pxc) • dS 



(6-9) 
(6-10) 



This relation is derived below by noting that the intrinsic internal energy is written 
as E = T s + E es + E xc . At a first step, we get the exchange-correlation pressure in 
the following way (see Appendix): 



SV 



dE x 



dV 



SN{n a } 



3V 
3V 



trP xc (ir 

dE xc 



dV 

dE xc 



SV 



5n(r) 



f SE XC 


dn(r) 


} 5n(r) 


SV dV 



(6-11) 
dr (6-12) 

SN{Tl a } 



dV 



Sn{r) 



+ / n(r)r • V 



5E X 



5n(r) 



dr . 



(6-13) 



SV 



with aids of (|4T1[) and (|5-7(l . Hence, on use of (|fi-13|) and (|5-8|) . the exchange- 
correlation pressure is derived as 



-SV 



dE x 



dV 



Sn(r) 



/ trP xc dr + / r-V-P xc dr 

Jv JV 



dv 



r-P xc -dS. (6-14) 



Since dE es /dV\ Sn ^ = and the kinetic pressure due to T s is given by (|4-6|) . the 
electron pressure is expressed finally as (j6TU|) with combined use of (|6T4j) . 

In the equilibrium, the occupation probability /, fulfills the following conditions 



dQ 

Wi 



TVnU(r) 



dF_ 

dfi 



TVNU(t) 



dE 

Wi 



SVNU(r) 



dfi 



0. 



(6-15) 



TVnoV ea {T) 



This can be verified below. In the expression of thermodynamic potential, Q 
T s — TS S + (-Ees + F xc ) ~~ A 4 -^, the derivative of the first term is written as 



dTs_ 

dfi 



dfi 



+ T — - + 

Tfi U eB d fi 



dn(r) 



with use of (|4-22j) . and the third term yields 



-Jf [Ees + ^xc] 



TVU(r) 



S[E es + F xc ] 



5n(r) 



TVU{r) 
[U cS (r) + ^}^dr . 



dn(r) 



dv 



(6-16) 

(6-17) 
(6-18) 
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because of d\E es -\- / dfi\ n (j.^ =0- Therefore, the derivative of Q with respect to fa 
becomes 



df2 

Wi 



TViiU(r) 



/\ . . An(v) dN 



TViM U eB (r) 



dN dn< 



TV M U Bf! (r) 



(6-19) 
, (6-20) 



which is equal to zero owing to (|4-21[) . Here, it is important to notice in this proof 
that the constant term in the effective potential U e g(r) should be — to fulfill (|6-15|) , 
Similarly, other relations in (|6-15|) can be proved with help of (|4-22|) . In this way, we 
can prove that the equilibrium conditions for the density distribution function n(r) 
and the Kohn-Sham eigenvalues to fulfill are given, respectively, by 



8Q 



5n(r) 



5F 



TV>l/(r) 



<5n(r) 



5E 



TVNU{r) 



and 



dQ 



7>[/(r) 



dF 

dei 



5n(r) 
dE 



5Q S 

SVNU(r) <M r ) 



0, (6-21) 



TVfM Q U cB (r) 



SNU(r) 



dei 



7> t/ cff (r) 



(6-22) 



TNU(r) 

Finally, we derive a finite-temperature generalization of the Janak theorem^ in 
the form: 

dE 

~~Tf~ = ti + Ml • (6-23) 

a Ji VU(r) 

which is different from original Janak's theorem in that the ground-state energy is 
replaced by the internal energy E in Eq. ()6-23j) . This can be easily proved with help 
of Eq. (|4-22|) , which is reduced to 



dfi 



v 



dS s f ,dn(r) 



because of df2 s /dfi\TVfj, u cft = 0. Since E = T s +E es +E XC: there results 



dE 



T 



VU(r) 



dS s 
dfi 



+ / 7eff(r 



dn(r) 



dr + 



6[E es + E xc ] 



Sn(r) 



SVU(r) 



dn(r) 



4 ln T^~T + A* = e i + Mi • 

P 1 - fi 

§7. Summary and discussion 



(6-24) 



dr (6-25) 



(6-26) 



In the first place, we described thermodynamics in the DF theory for a neutral 
electron- nucleus mixture in § [31 in a systematic manner. Next, we derived two funda- 
mental relations, (|4-22() and ()5-16|) . in the non-interacting system, which lead to the 
thermodynamical relations in the interacting system described in terms of quantities 
of the non-interacting system and the exchange-correlation effect. Thus, with com- 
bined of use of ()6-4j) , the relation (|4-22|) generates the pressure formula ()5-9|) as a force 
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theorem, the familiar Fermi distribution law^ (|6-15|) . a finite-temperature general- 
ization of Janak's theorem (|6-23j) and some other equilibrium conditions, (|6-21|) and 
Q6-22JI . which are translated into those of the non-interacting system, respectively. 
On the other hand, the relation Q5-16|l provides two types of force expressions, (|5-23|) 
and Q5-31J) (Andersen's force theorems ISIilSIJlQMa^ in conjunction with another type 
of force theorem (|6-1|) . which is equivalent to (|5-23|) . It is interesting to note that 
thermodynamic relations, (|2-15|) and (|2-16j) . provide two different but equivalent ex- 
pressions for the pressure: the one, (|6-8|) in the volume-integral form and the other, 
()6-lUj) in the surface-integral form, respectively. 

In the standard DF theoryP^'^nj^ the effective external potential is defined 

as 



5n(r) 



_ 5[E CS + E xc ] 
TVU(r) 6n ( r ) 



(7d) 

SVU(r) 



instead of the definition (|3-3U|) . where the interaction part of the chemical potential 
/xi is subtracted to give the effective external potential. Since the eigenvalues for 
£7 e fr(r) are related with those of U e s(r) as ej = ej + m, the formulas for f7 e fr(r) are 
obtained by replacing as ei—^iQ = ej— /i, and fio~ U e g(r) = fj,—U e g(r) from our results. 
For example, the occupation probability for a state Zi is given by /j = /(ej— /i) = 
/(ej— fJLo) with use of the chemical potential \x of the real system instead of hq. Also, 
the Janak theorem is written in the usual form: dE/dfi\y = e~j. On the other hand, 
the relation)^ e max = //, for the occupied maximum level e max at T = is expressed 
as e max = jUo for U cf£ (r). 

When we obtain the free-energy F of the electrons of a neutral electron-nucleus 
mixture on the basis of the DF theory with aids of the adiabatic approximation, the 
total pressure of this system is written aa^J 

3PV = 3(P e + P n )V 
dF 



- 3V 8V 



) +2(f n ) n + ^(R a -F a ) n (7-2) 

TN{R a }/ n 

= 2{f n ) n + 2(T s ) n + (E cs ) n + [ (trP xc ) n dr (7-3) 

Jv 

with the force on a-nucleus: F a = —V a F\TVN = —V a E es . Here, the bracket with a 
suffix n denotes the average concerning the nuclear variables defined as 

(A) _ Tr n {exp[-/?(f n + F)]A} 
,n Tr n {exp[-/3(f n + F)]} ' 

where Tr n means a trace using a complete set of eigenfunctions for the nuclear 
Hamiltonian, T n + FIRq,}, with the nuclear kinetic operator T n . In calculating the 
pressure for a solid, the nuclear pressure P n is usually neglected. 
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Appendix 



When F is a function of the independent variables {x\ , X2 , ■ ■ ■ , x n }, the derivative 
of F with respect to some parameter A is expressed as 



i=l,n 



dxi 
~dX 



(A-l) 



with Xi denoting a set of variables {x±, . . . , x n } except X{. Thus, there results 



dF 
lX 



X\ . ..Xm—\ 



dF 
IX 



+ 



dF 



dx r 



dx r . 



dX 



because of dF/dX\ xl ... Xjn = ^2i=m+i,n dF/dxi\ Xi ■ dxi/dX. 

On the other hand, when we transform the independent variables {x±, X2, 
to {xi,X2, ...,?/}, we get a formula akin to but different from l|A-2|) : 



OF 
dx\ 



__ d(F,x 2 , ...,x n ) / d(xi,x 2 , ...,y) _ 
Si d(xi,x 2 ,...,y) / d(xi,x 2 ,...,x n ) dxi 



+ 

xi.. .y 



dF 
dy 



dy_ 

Xn 



(A2) 
(A-3) 
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